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Preface 


A  laser  can  be  described  as  an  analog  computer.  For  example,  in  the  laboratory, 
the  inputs  to  a  laser  can  be  measured  as  a  voltage  and  current  applied  to  the  optical 
pump  source.  The  laser  operates  on  these  inputs  and  produces  an  output  that  can 
be  measured  as  an  energy  or  power.  By  observing  the  laser’s  output  for  various  input 
conditions,  a  mathematical  representation  of  the  laser’*  operation  can  be  developed. 
The  model  will  be  written  by  using  these  observations  iind  mathematical  conclusions 
to  represent  an  Nd:YAG/Glass  laser  on  a  digital  computer. 


Matthew  D.  Rotondaro 
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Abstract 

A  solid  state  laser  is  modeled  for  use  on  a  personal  computer.  The  model  can 
accurately  represent  Nd:YAG  and  Nd:Glass  lasers  during  both  continuous  wave  and 
repetitively  pulsed  operation.  The  large  number  of  input  variables  allow  the  user 
to  make  changes  to  the  cavity,  optical  pump  source,  and  the  cooling  system.  The 
model  calculates  the  energy,  power,  and  divergence  angle  of  the  beam  along  with 
other  variables.  Because  of  its  user-friendly  interface,  it  is  possible  for  someone  with 
limited  laser  experience  to  obtain  useful  information. 


CONTINUOUS  WAVE 
AND 

PULSED  SOLID  STATE  LASER  MODEL 


I.  Introduction 


1.1  Description  of  Model 

The  objective  of  this  thesis  is  to  model  a  solid  state  laser  during  both  repet¬ 
itively  pulsed  and  continuous  wave  operation,  for  both  stable  and  unstable  cavity 
designs.  There  were  three  important  considerations  during  the  construction  of  the 
model:  ensure  a  user-friendly  operating  environment,  model  the  laser  accurately, 
and  ensure  useful  and  easy  to  understand  outputs. 

The  model  is  broken  into  three  major  components:  laser  kinetics,  heat  transfer, 
and  wave  propagation.  Each  of  these  components  has  a  large  number  of  variable 
parameters,  giving  the  user  a  wide  variety  of  input  choices  (see  Appendix  B). 

The  model’s  most  important  capability  is  the  demonstration  of  transient  ther¬ 
mal  effects  of  a  repetitively  pulsed  laser.  For  example,  a  high-energy  pulsed  laser 
would  experience  changes  in  output  energy  per  pulse  as  a  function  of  pulse  rate  due 
to  heating  of  the  laser  rod.  The  model  could  be  used  to  make  predictions  about 
these  energy  changes.  In  addition,  the  model  could  also  examine  the  thermal  limits 
of  the  laser  rod  without  the  need  to  destroy  a  laser  rod. 

The  lasing  material  is  the  trivalent  rare  earth  Neodymium,  with  a  choice  of 
host  materials  including  YAG,  the  silicate  glass  ED-2,  and  the  phosphate  glass  Q-88. 
These  three  host  materials  are  well  documented  and  are  heavily  used  in  industry, 
making  them  good  choices  for  modeling. 
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The  software  was  written  using  Microsoft  QuickBASIC  because  of  its  effort¬ 
less  coding  environment.  Most  computers  have  some  implementation  of  the  BASIC 
language  making  the  code  portable  to  other  machines.  Also,  BASIC  is  such  a  well 
known  language  that  the  user  should  be  able  to  modify  the  code  to  suit  his/her 
personal  needs. 

1.2  Possible  Uses 

The  laser  model  has  a  wide  variety  of  possible  applications,  but  the  two  most 
interesting  in  terms  of  Air  Force  applications  are  laser  radar  and  blinding  optical 
sensors. 

Laser  radar  requires  a  high  power,  short  pulse,  single  frequency,  compact,  re¬ 
liable  laser.  The  compact  size  and  reliability  conditions  are  both  met  by  using  a 
diode  pumped  Nd:YAG  laser.  The  high  power  and  short  pulses  can  be  achieved  by 
Q-switching  the  laser.  The  laser  model,  with  its  ability  to  accurately  represent  a 
diode  pumped,  Q-switched,  repetitively  pulsed  Nd:YAG  laser,  could  be  quite  useful 
in  the  design  of  a  laser  radar. 

The  requirements  for  blinding  optical  sensors  are  very  similar  to  the  require¬ 
ments  for  a  laser  radar.  An  optical  sensor  can  be  temporarily  blinded  or  damaged 
by  the  application  of  a  high  power  pulse.  A  laser  used  to  blind  an  optical  sensor 
needs  a  high  pulse  repetition  rate  to  ensure  that  the  sensor  is  disabled  in  the  amount 
of  time  available  to  illuminate  the  target.  The  NdtYAG  laser  is  well  suited  to  this 
application,  but  it  only  produces  light  with  a  wavelength  around  l^m,  and  not  all 
sensors  are  in  this  region.  By  frequency  doubling  the  output  of  the  Nd:YAG  laser 
the  output  can  be  shifted  to  the  optical  band,  increasing  the  the  number  of  available 
targets  to  include  sensors  in  the  optica)  region.  The  tradeoffs  between  pulse  repeti¬ 
tion  rate  and  energy  deposited  per  pulse  on  the  target  can  be  evaluated  using  the 
laser  model. 
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II.  Background/Theory 


2.1  Four  Level  Laser 

Figure  1  shows  the  operation  of  an  optically  pumped  four-level  laser.  Initially, 
all  atoms  of  the  laser  material  are  in  the  lowest  level.  Excitation  is  provided  by 
optical  pumping  which  moves  ions  from  the  lowest  level  to  a  wide  absorption  band 
or  set  of  bands  labeled  level  three.  The  ions  populating  these  upper  bands  quickly 
depopulate  into  the  upper  laser  level,  level  two,  at  a  rate  of  where  r32  is  the 
relaxation  time  from  level  three  to  level  two.  The  laser  transition  or  stimulated 
emission  is  from  level  two  to  level  one.  Level  one  then  quickly  depopulates  down 
to  the  ground  level  at  a  rate  of  To  qualify  as  a  four-level  laser,  the  fluorescent 
lifetime  of  the  lower  laser  level  must  be  short  compared  with  the  fluorescent  lifetime 
of  the  upper  laser  level,  i.e.,  —  »  —.  In  addition,  the  lower  level  must  be  far 
enough  above  the  ground  level  to  prevent  it  from  becoming  thermally  populated. 

Figure  2  is  an  energy  level  diagram  for  neodymium  in  YAG  (5).  The  laser 
transition  with  the  best  performance  at  room  temperature  is  from  the  R?  line  to  the 
Y3  line,  making  i?2  and  Y3  the  upper  and  lower  laser  levels  and  4Fi  and  4Jn  the 
upper  and  lower  laser  manifolds.  The  population  of  a  level  inside  of  a  manifold  can 
be  described  using  the  Boltzmann  distribution: 

r  pie'w5 

A  ~ 

where  <?,-  is  the  degeneracy  of  level  i,  e,-  is  the  energy  of  level  i,  T  is  the  temperature, 
and  k  is  the  Boltzmann  constant. 

The  Boltzmann  factor  can  be  used  because  “•  •  •  It  is  also  assumed  that  the 
crystal  field  levels  within  each  manifold  are  in  quasi-thermal  equilibrium  at  all  times 
and  thus  their  populations  are  Boltzmann  distributed.”  (17).  For  example,  by  ob- 
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Figure  1.  Simplified  Energy  Level  Diagram  of  a  Four-Level  Laser 


taining  the  values  gx  =  2,  <72  =  2,  ej  =  0  cm-1  and  e 2  =  84  cm"1  from  Appendix  C, 
and  substituting  them  into  equation  1; 


2e~0.69J.6T 

2e°  +  2e~  o.ewor 

Then,  at  T  =  300  K,  /2  =  0:40.  This  means  that  40  percent  of  the  ions  populating 
the  level  two  manifold  can  be  found  in  the  i?2  level.  Note,  the  energy  of  the  lowest 
level  in  a  manifold  is  always  set  to  zero,  and  all  subsequent  levels  in  the  manifold 
are  measured  relative  to  the  lowest  level. 


2.2  Point  Model 

The  behavior  of  a  laser  can  be  represented  reasonably  well  by  using  a  set  of 
coupled  rate  equations  and  by  assuming  a  spatially  uniform  laser  medium.  In  the 
following  section,  the  rate  equations  will  be  laid  out  and  simplified  to  facilitate  a 
numerical  solution. 
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Figure  2.  Energy  Level  Diagram  for  Neodymium  in  YAG  (5) 


The  point  model  is  based  on  the  assumption  that  all  gains  and  losses  are  evenly 
distributed  throughout  the  laser  cavity.  More  specifically,  small  signal  gain,  favorable 
and  unfavorable  losses,  inversion,  and  photon  flux  are  all  constant  with  respect  to 
position  inside  the  cavity.  This  is  a  good  assumption  as  long  as  the  outcoupling  is 
not  so  large  as  to  create  a  substantial  gradient  in  gain  along  the  lasing  axis.  By 
using  these  assumptions  and  restrictions,  equations  3  to  7  can  be  used  to  describe  a 
four-level  laser  (17)  (5)  (7). 


dt 


n  gi  lc 


+  S 

Trc 


(3) 


dN3 

dt 


+  R3 


(4) 


dt 


-9—{N,  -  —Ni)fo  +  -  —  +  Rvh 

Tl  gi  TZ2  t2 


(5) 


dNi 

dt 


o-c.  g2  .  N2  f\  N3  iVi 

=  $—(N2 - Nx)fi  +  — —  +  — fi - +  Rifi 

n  gi  r2\  f2  r3 1  ti0 


(6) 


dN0 

dt 


Ni  N2 

Tiofl  T20/2 


T30 


—  R2  —  Rz 


(7) 


Equation  3  describes  the  change  in  the  number  of  photons  in  the  cavity  as  a 
function  of  time.  The  first  term  on  the  right  side  of  the  equation  is  the  increase  due  to 
stimulated  emission.  The  ratio  accounts  for  the  rod  being  shorter  than  the  cavity 
and  therefore  distributes  the  rate  of  stimulated  emission  across  the  whole  cavity. 
The  second  term  is  the  number  of  photons  lost  in  the  cavity  per  round  trip,  and  the 
last  term  is  the  rate  of  spontaneous  emission.  Note,  the  spontaneous  emission  term 
is  very  small  but  can  not  be  discarded  because  it  is  responsible  for  starting  the  lasing 
process. 
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Table  1.  Definition  of  the  Symbols  in  Rate  Equations  3  to  7 


$ 

Photons  per  unit  volume 

Nx 

Number  of  ions  in  level  x  where  x  =  0,1, 2, 3 

Rx 

Number  of  ions  excited  to  manifold  x  where  x  =  1,2,3 

9x 

Degeneracy  of  level  x  where  x  =  1,2 

u 

Boltzmann  factor  for  level  x  where  x  =  1,2 

B 

Rate  at  which  ions  depopulate  level  x  to  all 
lower  levels  where  x  =  1,2,3 

K293 

Rate  of  ion  depopulation  from  level  x  to  level  y 

B 

Time  it  takes  light  to  complete  one  round 
trip  inside  the  cavity 

lc 

Length  of  the  cavity 

mm 

Length  of  the  rod 

t 

Photon  losses  per  cavity  round  trip 

Stimulated  emission  cross-section 

n 

Index  of  refraction  of  the  laser  medium 

Rate  of  spontaneous  emission  into  the  0,0  transverse  mode 

n 

Speed  of  light  in  a  vacuum 

Equation  4  describes  the  change  in  the  number  of  ions  excited  to  the  level  three 
manifold  as  a  function  of  time.  The  first  term  on  the  right  side  of  the  equation  is 
the  loss  from  the  level  three  manifold  to  all  lower  levels,  through  phonon  relaxation. 
The  second  term  is  the  increase  due  to  absorption  of  energy  by  the  ions  from  the 
optical  pump  source. 

Equation  5  describes  the  change  in  the  number  of  ions  populating  the  upper 
laser  level  as  a  function  of  time.  The  first  term  on  the  right  side  of  the  equation  is  the 
loss  from  level  two  due  to  stimulated  emission.  The  second  term  is  the  increase  to 
level  two  from  the  pump  bands  (level  three)  through  phonon  relaxation.  The  third 
term  is  the  loss  from  level  two  to  all  lower  levels  through  spontaneous  emission  and 
through  phonon  relaxation.  The  last  term  is  the  increase  due  to  the  optical  pump 
source.  The  Boltzmann  factor  /2  is  necessary  because  a  change  of  one  ion  in  the 
level  two  manifold  only  changes  the  upper  laser  level  by  the  amount  fi- 

Equation  6  describes  the  change  in  the  number  of  ions  occupying  the  lower 
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laser  level  as  a  function  of  time.  The  first  term  on  the  right  side  of  the  equation 
is  the  increase  in  the  lower  level  due  to  stimulated  emission.  The  second  and  third 
terms  are  the  increases  to  level  one  from  levels  two  and  three  through  spontaneous 
emission  and  phonon  relaxation.  The  Boltzmann  factor  /2  is  in  the  denominator  of 
the  second  term  to  account  for  the  increase  from  manifold  two,  not  just  the  increase 
from  level  two.  The  fourth  term  is  the  loss  from  level  one  to  the  ground  state  due  to 
phonon  relaxation  and  the  last  term  is  the  increase  due  to  the  optical  pump  source. 
The  Boltzmann  factor  fi  is  necessary  because  a  change  of  one  ion  in  the  level  one 
manifold  only  changes  the  level  by  the  amount  fi. 

Equation  7  describes  the  change  in  the  number  of  ions  populating  the  ground 
level  as  a  function  of  time.  The  first  three  terms  are  the  ions  that  eventually  make 
it  back  down  to  the  ground  level  after  being  excited  to  the  upper  manifolds.  The 
last  three  terms  are  the  losses  to  the  ground  state  level  due  to  optical  pumping  of 
the  upper  manifolds.  The  Boltzmann  factors  f\  and  /2  are  included  in  the  first  and 
second  terms  to  account  for  the  increase  to  the  ground  level  from  manifolds  1  and 
2,  not  just  from  levels  1  and  2. 

These  five  rate  equations  are  really  more  than  is  needed  to  explain  what  is 
going  on  during  the  laser  operation.  The  first  simplification  is  the  elimination  of 
equation  7.  For  Nd:YAG  and  Nd:Glass  lasers  only  a  small  fraction  of  the  total 
number  of  ions  in  the  ground  state  are  ever  moved  to  the  higher  levels-usually  on 
the  order  of  only  1  percent.  This  implies  that,  during  the  entire  lasing  process,  the 
total  number  of  ions  in  the  ground  state  is  approximately  equal  to  the  total  number 
of  ions.  Because  the  number  of  ions  in  the  ground  state  does  not  change  significantly 
then  the  absorption  of  pump  energy  from  the  ground  state  ions  remains  constant 
throughout  the  lasing  process  and  the  number  of  ions  in  the  ground  state  does  not 
need  to  be  tracked  throughout  the  lasing  process. 

The  next  simplification  is  a  bit  more  rigorous.  Equation  4  can  be  neglected 
when  t32  <  PumP  Ration .  gqUatjon  4  has  a  solution  of  the  form 
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(8) 


N3  =  #373(1  -  exp(— )) 
73 


where 


_1_ 

7-3 


730 


1 

- h 

731 


_1_ 

732 


(9) 


and 


732  731 ,  T30 


Therefore 


(10) 


N3  «  #3X32(1  -  exp(— ))  (11) 

732 

This  equation  shows  that  when  t  is  greater  than  1  /is  and  T32  is  100  ns  or 
smaller  then 


N3  ~  #3X32  . 


(12) 


where  t  is  the  amount  of  time  the  flash  lamp  is  on  and  is  typically  greater  than  1  /is. 
Equation  12  can  be  substituted  directly  into  equation  5  to  yield 


^  (JV,  -  +  Rzh  -  —  +  Rzh 

at  n  gi  t2 

and  by  letting  #2  =  #2  +  #3  then 


(13) 


dN2  ,<rc  02,rw  N2  D,  , 

=  -$—{N2  -  —Ni)f2 -  —  +  R2h 

at  n  <7i  X2 


(14) 
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This  means  that  if  r32  is  100ns  or  smaller,  then  equation  4  can  be  eliminated 
and  the  manifold  pump  rate  R3  can  be  directly  added  to  the  manifold  two  pump 
rate. 

The  pump  rates  are,  in  general,  very  difficult  to  calculate.  They  require  the 
modeling  of  both  the  optical  pump  source  and  the  absorption  of  the  lasing  material. 
Such  an  undertaking  would  be  a  considerable  amount  of  work  and  would  be  outside 
the  scope  of  this  project.  Therefore,  the  following  approximation  is  used  to  calculate 
the  pump  rate  R'2. 


g  _  Electrical  power  into  the  lamp  xThe  fraction  for  lasing  . 

2  h  x  Laser  frequency  x  Volume  of  the  rod  ' 

where  Electrical  power  into  the  lamp  and  The  fraction  for  lasing  are  user  inputs  (see 
Appendix  B)  and  h  is  Planck’s  constant.  This  is  a  good  approximation  for  the  Nd 
ion  lasers  because,  as  stated  earlier,  the  absorption  cross-section  can  be  treated  as 
a  constant.  Therefore,  the  pump  rate  is  constant  and  determined  only  by  a  fraction 
of  the  flash  lamp  input  electrical  power. 

The  model  is  designed  to  accurately  represent  the  output  of  a  pulsed  laser. 
A  subset  of  the  pulsed  laser  is  the  Q-switched  laser.  The  Q-switch  is  a  device 
placed  inside  the  laser  cavity  and  is  generally  used  to  generate  high  pulse  powers. 
The  Q-switch  prevents  the  laser  from  lasing  by  spoiling  the  quality  of  the  cavity. 
This  is  done  so  that  the  population  inversion  can  build  past  the  minimum  inversion 
required  to  lase;  the  critical  inversion  Ncrit-  Finally,  by  the  time  the  Q-switch  is 
opened  and  lasing  is  allowed  to  take  place,  the  inversion  has  increased  well  above  the 
critical  inversion.  This  large  inversion  produces  a  large  gain  in  the  cavity  causing 
the  number  of  photons  to  build  very  rapidly.  The  rapid  buildup  in  photons  extracts 
the  energy  from  the  inversion  very  quickly  producing  the  desired  high  power  output 
pulse. 
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While  the  quality  of  the  cavity  is  spoiled  by  the  Q-switch,  the  number  of  pho¬ 
tons  in  the  cavity  is  approximately  equal  to  zero.  Therefore  $  =  0  and  equation  14 
can  be  written  as 


dN2 

dt 


+  B!2h 


which  has  the  following  solutions  for  the  specified  domains 


(16) 


N2  =  ^72/2(1  -  exp(— ))  0  <  t  <  tQ-axuitch  <  T  (17) 

t2 

and 

N2  =  N2(T)  exp(— )  T  <  t  <  tQ.awitch  (18) 

t2 

where  T  is  the  time  at  which  the  flash  lamp  was  turned  off,  and  tQ~awitch  is  the  time 
that  Q-switching  occurred. 

Even  after  these  simplifications,  equations  3  to  7  are  still  a  set  of  nonlinear 
coupled  differential  equations  which  require  n  .merical  techniques  to  solve.  In  order 
to  apply  numerical  techniques,  it  is  advisable  to  first  normalize  the  equations. 

A  quantity  well  suited  to  the  normalization  of  the  rate  equations  is  the  critical 
inversion.  The  critical  inversion  is  the  minimum  inversion  necessary  to  achieve  lasing 
and  is  defined  during  steady  state  operation.  Steady  state  means  that  the  net  change 
in  the  number  of  photons  per  unit  time  in  the  cavity  is  equal  to  zero: 


n 

dt~° 


Equation  3  may  now  be  written  as 


—  =  — (jV2  -  +  S 

Trc  n  Pi  lc 


(19) 
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The  spontaneous  emission  rate  S  is  given  by  (7) 


N2  c3 
Til  87TI/2Ai 'Vcavily 


(21) 


where 

v  =  The  laser  frequency, 

A v  =  FWHM  of  the  lineshape  for  atomic  response, 

Vcavity  =  volume  of  the  laser  cavity. 

Noting  that  S  is  very  small  and  may  be  neglected  allows  the  critical  inversion  to  be 
defined  as 


AT0ril  =  (Ni  -  — JVi)  = 

Si 


en  lc 
rTCac  lr 


(22) 


Another  parameter  useful  for  normalizing  the  rate  equations  is  the  cavity  life¬ 
time  (7). 


Cavity  round  trip  time 

cav  —  J?rac^,ono|  round  -trip  loss 

The  cavity  life  time  for  a  stable  cavity  is  given  by  (7) 


(23) 


Tcav  —  ~  [fib  (^c  ^r)] 

c 


ln(flfe)  +  2a/r  +  ln(J) 


Q'J 


(24) 


Table  2.  Definition  of  Terms  for  the  Stable  Cavity  Lifetime 


Ri 

Reflectivity  of  mirror  one  =  1  -  fractional  outcoupling 

R* 

Reflectivity  of  mirror  two  =  1  -  absorption  losses 

Losses  at  the  mirrors 

2  air 

Losses  due  to  the  rod 

Q 

Transmission  of  the  Q-switch 

Ms) 

Losses  due  to  the  addition  of  a  Q-switch  into  the  cavity 

The  cavity  life  time  for  an  unstable  cavity  is  given  by 
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1 


(25) 


Tcov  —  —  [nlr  +  (lc  ~  L)] 
C 


^ r  +  ^(q)] 


Table  3.  Definition  of  Terms  for  the  Unstable  Cavity  Lifetime 


Rx 

Reflectivity  of  mirror  one  =  1  -  absorption  losses 

Ri 

Reflectivity  of  mirror  two  =  1  -  absorption  losses 

^n(«ikr) 

Losses  at  the  mirrors  and  unstable  cavity  outcoupling 

2  alr 

Losses  due  to  the  rod 

Q 

Transmission  of  the  Q-switch 

r 

1  -  Unstable  cavity  outcoupling 

Mb) 

Losses  due  to  the  addition  of  a  Q-switch  into  the  cavity 

e  = 


(26) 


Ncrit  — 


n  lc 
(TCTCav  lr 


(27) 


Now,  normalization  can  proceed  by  using  the  following  expressions: 


Merit 


nx 


Nx 

Ncrit 


t 


x,y  — 


Tcav 


=  Rx 


Merit 


(28) 

(29) 

(30) 

(31) 

(32) 
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(33) 


s  =  S 


Ncrit 


Finally  writing  the  normalized  equations  as 


^  =  <}>{n2  -  — »i)  -<f>  +  s 
at  gi 


(34) 


dn3  n3  , 

ir  =  ~i;+r 3 


(35) 


dl^2  u  52  tc  ,  n3  t  n2  .  r 

—  =  —<j>(n2 - nO/27-  +  7-/2  -  —  +  ^2/2 

at  g  1  /r  i32  12 


(36) 


d«i  /  /  52  w  4  .  «2  /1  n3  »1  , 

—  =  ^(n2  -  — ni)/i—  +  — y  +  — /i  -  — •  +  ri/i 

at  g  1  <r  121  /2  t3i  tio 

and  for  the  case  when  r32  is  of  the  order  of  100ns  or  smaller,  then 


(37) 


$  =  ^("j  -  — n,)  -  ij>  +  s 
at  g  1 


(38) 


dn2  , ,  52  w  *c  TI2  ,  . 

__  =  _^(n2__ni)/2___+r2/2 


(39) 


92n  w  /c  n2  /1  r3<32  .  ni 

— -  =  <p(n2  -  --nijfiy  +  — -r  +  - — /1  ~  +  n  /1 

at  5l  *r  t21  J2  t3i  IlO 


(40) 


The  last  thing  that  needs  to  be  determined  before  integration  can  proceed  is 
the  stimulated  emission  cross-section  a.  The  stimulated  emission  cross-scction  is 
given  by  (5)  as 


a  =  A21 


A2 

87m2 


M 


(41) 
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where 

An  —  The  Einstein  A-coefficient 
A  =  The  laser  wavelength 
n  =  The  index  of  refraction  in  the  rod 
g(v)  =  The  atomic  line  shape 

Thermal  broadening  is  the  mechanism  responsible  for  the  linewidth  of  the 
Nd:YAG  laser.  “Thermal  broadening  is  brought  about  by  the  effect  of  thermal 
lattice  vibrations  on  the  atomic  transition.  The  thermal  vibrations  of  the  lattice 
surrounding  the  active  ions  modulate  the  resonance  frequency  of  each  atom  at  a  very 
high  frequency.  This  frequency  modulation  represents  a  coupling  mechanism  between 
the  atoms,  therefore  a  homogeneous  linewidth  is  obtained.”  (7)  The  homogeneous 
linewidth  leads  to  a  Lorentzian  line  shape  for  atomic  response,  such  that 


9 (")  =  (^)  (*  “  *0’  +  ( 


Av 


\2 


-1 


(42) 


The  laser  model  assumes  lasing  only  at  the  central  peak  v0  and  therefore  the  line 
shape  can  be  reduced  to 


M  = 


2 

7tAi/ 


(43) 


An  inhomogeneously  broadened  line  shape  occurs  in  the  fluorescence  of  a 
neodymium  doped  glass.  “There  are  variations  from  rare  earth  site  to  rare  earth 
site  in  the  relative  atomic  positions  occupied  by  the  surrounding  lattice  ions.  This 
gives  rise  to  a  random  distribution  of  static  crystalline  fields  acting  on  the  rare  earth 
ions.  Since  the  line  shifts  corresponding  to  such  crystal-field  variations  are  larger, 
generally  speaking,  than  the  width  contributed  by  other  factors  associated  with  the 
transition,  an  inhomogeneous  line  results. ”(7)  The  inhomogeneous  line  can  be  rep¬ 
resented  by  a  Gaussian  frequency  distribution  given  by 
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Again,  the  laser  model  assumes  lasing  only  at  the  central  peak  Uo  and  therefore  the 
line  shape  can  be  reduced  to 


Now  equations  38  to  40  can  be  solved  by  using  a  Runga-Kutta  or  an  Adams 
Predictor- Corrector  method  of  integration  (see  Appendix  A).  The  integration  is 
needed  to  determine  the  number  of  photons  in  the  cavity  at  any  given  time.  Once 
the  number  of  photons  in  the  cavity  is  known,  then  the  number  of  photons  in  the 
beam  can  be  determined  by  first  rewriting  rcau  in  equation  3  in  terms  of  favorable 
and  unfavorable  losses. 

*££W -**,)£:  _JL__*_  +  S 

dt  7i  g\  lc  Tjav  Tunjav 

Favorable  losses  are  the  losses  that  go  into  the  laser  output,  and  unfavorable  losses 
include  losses  in  the  rod,  Q-switch,  and  the  back  mirror. 


'T’cav  Tfav  Tun/av 


Separating  out  the  favorable  losses  from  equation  46, 


(47) 


d^beam  ^ 

dt  Tcav 

This  is  the  number  of  photons  in  the  beam  per  unit  volume  per  unit  time,  where 
Tjav  for  a  stable  cavity  is  given  by 
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Figure  3.  Geometry  of  the  Laser  Rod 


Tfav  —  [u/r  -j-  (/c  /r)] 

c 

and  Tjav  for  an  unstable  cavity  is  given  by 


(49) 


T/au  =  ~  [n<r  +  (4  -  M] 
c 


(50) 


2.3  Nodal  Network  Heat  Transfer 

An  inherent  effect  of  shining  a  light  on  the  laser  rod  to  pump  the  laser  in¬ 
version  is  the  heating  of  the  rod.  This  heating  is  responsible  for  thermal  lensing, 
birefringence,  stress  fractures,  and  in  extreme  cases,  melting  the  rod.  Therefore, 
a  description  of  the  temperature  distribution  in  the  rod  is  necessary  to  accurately 
describe  the  laser.  In  this  section  the  temperature  distribution  in  the  rod  will  be 
found  by  solving  equation  51  which  describes  the  heat  flow  in  a  flash  lamp  pumped, 
fluid  cooled  rod  (2). 


<h,(T)p(T)-^  =  V  •  (k(T)VT)  +  r,  0,  z)  (51) 

The  temperature  T  =  T(r,  0 ,  z ,  t)  and  is  subject  to  the  following  boundary  conditions: 
T(O,0,  z,t)  <  oo 
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r(r,  0, 2,  t)  =  T(0, 27 r,  2,  t) 

=  *f*(r(B, »,*,«)  -  T, i,„) 
^),=o  =  *f-(T(r,e,0,t)  -  Tait) 

=  k-t(T(r,6,L,t)  -  Tair) 
T(r,0,2,  0)  =  r//ut(i 


Table  4.  Definition  of  Terms  for  the  Heat  Equation 


T 

Temperature  distribution  in  the  rod 

cp 

Heat  capacity  of  the  rod 

P 

Density  of  the  rod 

k 

Thermal  conductance  of  the  rod 

9* 

Heat  energy  entering  the  rod  per  unit  volume 

A 

Surface  area  normal  to  r .  See  figure  3 

B 

Surface  area  normal  to  z.  See  figure  3 

TjUm 

Temperature  of  the  fluid  film  surrounding  area  A 

Tair 

Temperature  of  the  air  surrounding  area  B 

Equation  51  is  very  difficult  to  solve,  and  is  more  general  than  is  required  to 
do  an  accurate  job  of  describing  the  heat  flow  in  the  laser  rod.  The  equation  can 
be  simplified  by  eliminating  the  9  and  z  dependence.  This  is  accomplished  by  using 
the  following  three  assumptions:  q*  =  Constant ,  the  heat  flow  at  the  ends  of  the 
rod  is  small,  and  the  fluid  flow  on  surface  A  is  fast  enough  to  ensure  that  the  fluid 
temperature  does  not  change  significantly  as  it  travels  the  length  of  the  laser  rod. 
Using  these  assumptions,  equation  51  can  be  reduced  to 

<h{T)p(T)^  =  V  •  (k(T)VT)  +  q-  (52) 

.The  temperature  T  =  T(r,t)  and  is  subject  to  the  following  boundary  conditions: 
T(0,t)  <  00 

T(r,0)  =  T fiM 

This  can  be  further  simplified  by  assuming  that  cp,  p  and  k  are  constants,  yielding 
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Figure  4.  One  Dimensional  Heat  Transfer  Nodal  Network 


|  =  Av2r+-^-  '  (53) 

dt  Cp  p  Cp  p 

The  temperature  T  —  T(r ,  t)  is  subject  to  the  following  boundary  conditions: 

T(0,t)  <  oo 

T(r,  0)  =  T fluid 

The  heat  q *  is  given  by 

„  _  Flash  Lamp  Wattage  x  Power  Fraction  Heat  .  . 

^  Laser  Rod  Volume 

where  Flash  Lamp  Wattage,  Power  Fraction  Heat,  and  Laser  Rod  Volume  are 
all  user  inputs  (see  Appendix  B). 
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An  approximate  solution  may  be  had  by  using  a  finite  difference  approximation 
(see  Figure  4)  (2).  The  solution  to  equation  53  can  be  written  as 


T-  —  T{  6t 


y'  Ti _ 1i  4.  Hi. 

y  KijCi  *  ci 


(55) 


where  the  thermal  resistance  between  nodes  i  and  j  for  conduction  inside  the  rod  is 


_  Sr 
nii=kA 


(56) 


'J 


the  thermal  resistance  between  nodes  i  and  j  for  convection  at  the  rods  surface  is 


(57) 


and  the  thermal  capaciy  for  node  i  is  given  by 


C{  —  Vipcp 


(58) 


Table  5.  Definition  of  Terms  for  the  Finite  Difference  Solution  to  the  Heat  Equation 


Tt 

Temperature  at  node  i 

T[ 

Temperature  at  node  i  at  a  time  t  +  St 

Heat  input  at  node  i 

Sr 

Distance  between  nodes 

St 

Time  step 

h 

Heat  transfer  coefficient  for  the  cooling  fluid 

A-jj 

Surface  area  evaluated  half-way  between  node  i  and  j 

Ci 

Thermal  capacity  for  node  i 

Vi 

Volume  of  node  i 

P 

Density  of  the  rod  material 

Cp 

Heat  capacity  of  the  rod  material 

“These  equations  just  given  express  the  temperature  at  a  given  node  in  terms 
of  a  discrete  spatial  step  Sr  and  discrete  time  step  St.  The  important  point  is  that 
the  continuous  solution  of  the  original  differential  equation  is  replaced  by  a  numerical 
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Temperature  K 

Figure  5.  Experimental  vs  Cubic  Curve  Fit  for  Air 


solution  of  these  difference  equations  at  particular  points  in  space  and  at  particular 
intervals  of  time.”  (2) 

The  solution  is  only  stable  if  8t  is  selected  such  that  (2) 

St  < - \ -  (59) 

The  smallest  8t  for  all  nodes  is  the  one  necessary  to  insure  stability  of  the  solution 
(2).  “The  implication  of  equation  59  is  that  the  choice  of  the  time  step  8t  is  inti¬ 
mately  connected  with  the  choice  of  the  spatial  increment,  6r ,  which  is  involved  in 
the  resistance  7 Zij.  As  smaller  spatial  increments  are  chosen  to  reduce  the  associ¬ 
ated  truncation  error,  smaller  time  increments  must  also  be  used  in  order  to  satisfy 
equation  59.  Thus,  increased  accuracy  in  the  spatial  network  is  obtained  at  the  cost 
of  smaller  time  increments”  (2). 

By  using  a  one-dimensional  nodal  network  (see  Figure  4),  the  radial  temper¬ 
ature  distribution  in  the  rod  can  be  found.  Initially,  the  temperature  at  all  of  the 
nodes  is  equal  to  the  temperature  of  the  water  Tw  surrounding  the  laser  rod.  When 
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Figure  6.  Experimental  vs  Cubic  Curve  Fit  for  Water 


the  rod  is  heated  through  flash  lamp  pumping,  it  is  assumed  that  the  rod  heats  up 
uniformly.  As  the  temperature  of  the  laser  rod  increases  a  temperature  gradient  be¬ 
tween  the  rod  and  the  water  develops  and  heat  begins  to  flow  out  of  the  rod  creating 
a  parabolic  temperature  distribution  in  the  rod. 

The  nodal  network  solution  allows  for  generalization  to  the  more  complex  cases 
where  the  heat  capacity  and  thermal  conductivity  are  not  constants.  Another  ad¬ 
vantage  to  reducing  the  more  general  case  with  the  variable  conductivity  and  heat 
capacity  is  that  the  steady  state  solution  is  easily  evaluated  and  can  be  used  to 
validate  the  numerical  solution  before  generalizing  to  the  more  complex  case. 

At  the  surface  of  the  rod  the  thermal  resistance  is  given  by  equation  57.  This 
resistance  is  dependent  on  the  heat  transfer  coefficient  h.  The  heat  transfer  coefficient 
is  a  function  of  temperature  and  velocity.  Figures  5  and  6  are  plots  of  hyjfy  verses 
temperature  for  dry  air  and  water  (2)(6).  The  second  curve  in  each  figure  is  a  third 
degree  polynomial  curve  fit. 

hJ pr  =  4.04562-5.82881o:10-4T/,7m+1.92359a:10-7T/2t-,m-1.98457xl0-11T/3j7m  (60) 
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for  air,  and 


=  -3058.2  +  23.89lT/,/m  -  .0453037^  +  2.80149xl0-5T/3i,m  (61) 

for  water.  L  is  the  length  of  the  laser  rod,  and  V  is  the  velocity  of  the  coolant  flowing 
over  the  rod. 

When  a  fluid  of  one  temperature  flows  over  a  rod  at  a  different  tempera¬ 
ture,  a  thermal  boundary  layer  develops.  Good  results  are  obtained  when  the 
fluid  properties  are  evaluated  at  the  film  temperature  Tjum  given  by  the  following 
approximation  (6) : 


T/ilm  = 


Tjluid  +  TWall 


(62) 


The  finite  difference  solution  is  general  enough  to  describe  the  temperature 
distribution  in  the  rod  at  any  time  t.  The  problem  with  the  finite  difference  solution 
is  that  as  t  — ►  oo,  the  computational  time  also  goes  to  infinity.  This  condition  is 
encountered  by  a  continuous  wave  laser.  The  resolution  to  this  problem  can  be  found 
by  recognizing  that  equation  53  has  an  analytic  solution  in  the  steady  state  limit. 
See  Chapter  3  for  details. 

The  problem  with  the  analytic  solution  is  that  it  depends  on  the  heat  transfer 
coefficient  which  is  a  function  of  temperature.  Because  of  this  dependence,  the 
analytic  steady  state  solution  is  really  a  transcendental  equation.  Using  the  steady 
state  solution  from  Chapter  3, 


T(R)  = 


T fluid  +  T(R) 
2 


q*  2  kR 

+  ik  {h(T(R)) 


(63) 


By  solving  this  transcendental  equation  for  T(R),  and  by  using  the  steady  state 
solution  in  Chapter  3,  the  steady  state  temperature  distribution  in  the  rod  can  be 
found. 
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2.4  Transient  Thermal  Effects 

As  the  rod  changes  temperature  the  index  of  refraction  also  changes.  This 
change  in  index  can  be  described  by  the  following  expression: 

dn 

«.■=»-  POO  -  T( 0)]—  (64) 

where  n0  is  the  index  of  refraction  at  the  center  of  the  rod,  n  is  the  index  of  refraction 
of  the  material  at  300°if,  T(0)  is  the  temperature  at  the  center  of  the  rod,  and  is 
the  change  in  index  of  refraction  with  respect  to  temperature.  The  index  of  refraction 
as  a  function  of  radius  is  given  by 

n(r)  =  no+[T(r)-T(0)]^  (65) 

This  change  in  index  accounts  for  the  thermal  lensing  in  the  rod.  Note,  for 
example,  that  if  ~j~  is  positive,  and  the  rod  is  hotter  at  the  center  than  on  the  walls, 
then  the  rod  will  act  like  a  positive  lens. 

2.5  Beam  Quality 

The  model  can  be  used  to  analyze  a  laser  with  a  stable  or  an  unstable  cavity. 
Knowing  the  cavity  design  is  very  important  because  the  techniques  for  determining 
information  about  the  beam  are  very  different  for  the  two  cases.  In  this  section  the 
method  of  determining  the  stability  of  a  cavity  will  be  derived.  This  will  be  followed 
by  a  discussion  of  the  techniques  used  for  determining  spot  size,  mode  volume,  and 
divergence  angle  for  a  stable  cavity.  The  section  will  then  conclude  with  a  discussion 
of  the  techniques  used  to  determine  the  outcoupling,  spot  size,  and  divergence  angle 
for  an  unstable  cavity. 

2.5.1  Stability  A  cavity  is  considered  stable  if  a  ray  inside  the  cavity  can 
go  through  an  infinite  number  of  round  trips  and  still  remain  inside  the  cavity.  If 
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Figure  7.  The  Cavity  Design  Used  by  the  Laser  Model 


Unit  Cell 


Figure  8.  Equivalent  Cavity  Lens  Wave  Guide 
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the  rays  were  able  to  ’’walk”  out  of  the  cavity  then  the  cavity  would  be  considered 
unstable.  Stability  can  be  found  by  using  transformation  matrices  that  transform  an 
input  ray  at  position  z  to  the  position  z’  (see  Figure  9).  The  transformation  matrices 
are  constructed  based  on  a  paraxial  approximation;  tan(0)  «  sin(0)  «  6  =  r'. 
Figure  7  depicts  the  cavity  design  used  by  the  laser  model.  Figure  8  is  the  equivalent 
lens  wave  guide.  The  following  matrices  describe  the  round  trip  propagation  of  a  ray 
through  the  cavity  (12)  (9)  (8)  (5). 


rT 

'  1  let'1 

1  °) 

l 

.  rT  . 

1  / 

vfe:  n) 

\ 

cos(6)  ^ b  sin(6) 

~zSsin(6)  cos(&) 


1  0 

n— 1  1 


nRrc 

1  0 
1— n 


\ 

(i 

Let ' 

'  1  o\  ' 

/ 

1  / 

\  Rm  1  /  V 

\  Rr 


n 


^sin(&) 


cos  (b)  L 

{  -Z^Mb)  cos (6)  ) 


\  nR, 


where 

n  =  Index  of  refraction  averaged  across  the  radial  direction, 
no  =  Index  of  refraction  at  the  center  of  the  rod. 
hr  =  Index  of  refraction  at  the  surface  of  the  rod. 


(66) 
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and  where 


7 ]R2 

nR  ~  n0 - — 

(67) 

2n0  +  nR 
"  =  3 

(68) 

2(n0  -  nfl) 

''  B? 

(69) 

V  no 

(70) 

cos(6)  i$Mb)y 

f^sin  (b)  cos  (6)  j 

(71) 

is  the  transformation  matrix  for  a  laser  rod  acting  like  a  positive  lens. 

When  hr  >  no,  the  rod  acts  like  a  diverging  lens,  and  b  becomes  imaginary.  By 
letting  b  =  iB  where  B  =  LrodyJ —jfc  and  using  some  very  well  known  trigonometric 
identities  [ Cos(iO )  =  Cosh(9)\  and  [ Sin(i9 )  =  iSinh(9 )]  the  transformation  matrix 
for  a  laser  rod  acting  like  a  diverging  lens  can  be  written  as 


c°sh  (B)  ^sinh(B)| 

V  &inh(£)  cosh(B)  ) 

By  multiplying  the  matrices  in  equation  67  together  from  right  to  left,  equa¬ 
tion  67  can  be  written  as 


.4  B 
C  D 
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Completing  the  multiplication,  the  following  recurrence  relations  can  be  written: 


rr+ i  =  Arp  +  Br'T  =>  r'T  =  -^(rr+i  —  Arp) 

(74) 

rr+i  =  Ctt  +  Dr'T 

(75) 

rT+ 1  =  ^(rr+2  ~  Arp+i)  =  Crp  +  Drip 

(76) 

— (rr+2  -  Arp+i)  =  Crp  +  —  (rr+i  -  Arp) 

(77) 

rr+ 2  2  2  rp+i  +  (AD  BC)rp  =  0 

(78) 

This  second-order  difference  equation  describes  the  rays  that  pass  through  any 
cavity  represented  by  an  ABCD  matrix.  The  solution  will  be  found  by  assuming  a 
solution  of  the  form  rp  =  r0et6T.  Substituting  this  solution  into  equation  78  yields 


r0e 


i$T 


J28 


A  -f  D 


-  +  (AD  -  BC)\  =  0 


'J 


(79) 


In  order  for  there  to  be  a  non-trivial  solution,  the  expression  inside  of  the  brackets 
must  be  equal  to  zero.  Therefore,  e'9  can  be  found  by  using  the  quadratic  equation: 


e„  =  A  +  D±. 


AD-BC- 


(A  +  D' 

V  2  , 


(80) 


In  order  for  the  solution  to  be  oscillatory, 

AD-BC- 


(81) 


and  therefore,  to  ensure  stability  inside  the  cavity: 


\A  +  D\  <  2 VAD  -  BC 


(82) 
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2.5.2  Stable  Cavity  Once  the  stability  of  the  cavity  has  been  determined, 
the  spot  size  of  the  beam  everywhere  along  the  lasing  axis  must  be  found.  This  will 
allow  the  determination  of  the  mode  volume,  and  the  divergence  angle  of  the  beam. 
In  order  to  find  the  spot  size  for  a  stable  cavity  it  will  be  necessary  to  introduce  the 
complex  beam  parameter  q  (5): 

I- I_ •  A 

q  R  miw2 

where 

R  =  radius  of  curvature  of  the  wave  front 
w  —  radius  of  the  spot  defined,  as  the  point  where  the  intensity  falls  off  to  of  the 
peak  intensity  for  the  0,0  transverse  mode. 

If  the  radius  R  and  the  spot  size  w  are  known  at  some  position  inside  the 
cavity,  then  it  is  possible  to  find  the  complex  beam  parameter  at  any  other  point 
inside  the  cavity  by  using  the  expression 

(84) 

where  the  matrix 


I  =  £±i 

92  A+  f 


A  B 
C  D 

is  a  transformation  matrix  from  the  point  of  known  q  to  the  point  where  q  needs 
to  be  determined.  At  a  mirror  of  a  stable  cavity,  the  radius  of  curvature  of  the 
wavefront  is  the  same  as  the  radius  of  curvature  of  the  mirror.  Also,  for  a  round  trip 
in  a  stable  cavity,  q\  =  q^.  Therefore  defining 


(86) 
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and 


_  AC  +  (AD  +  BC)xx  +  BD{x\  +  yy) 
X 2  A2  +  2ABxi  +  B2(xl  +  yl) 

yi(AD  -  BC) 

V 2  A2  +  2  ABxx  +  B2{x\  +  y\) 
at  a  mirror  where  xi  =  x2  =  x  and  yx  =  y^  =  y,  then 

\AD  -  BC  -A2-  2ABx  .1> 
«=[ - g5 - 1 


(89) 

(90) 

(91) 


Now  that  q  is  known  at  a  mirror,  then  it  can  be  propagated  through  the  cavity 
to  obtain  the  spot  size  throughout  the  cavity.  If  the  spot  is  propagated  through 
the  rod  using  a  small  step  size,  Sz,  then  the  spot  size  everywhere  in  the  rod  can  be 
determined.  The  mode  volume  can  be  calculated  by 

N 

W  =  f>zJ2  TTWi  (92) 

1=1 

where  6z  =  R?iL™ath ,  The  mode  volume  is  for  the  (0,0)  cylindrical  geometry  trans¬ 
verse  mode,  Vo.o  =  ^Modt-  The  spot  size  for  higher  order  transverse  modes  can  be 
found  by  using  Wpj  =  CPi/tt>o,o  (7),  and  therefore  the  higher  order  mode  volume 
can  be  found  using  VPii  =  C^Vo, o,  where  CPii  is  a  constant  for  the  mode  p,l.  For 
example,  Co,o  =  1,  Con  =  1.5,  Ci,o  =  1.9,  and  C\,\  =  2.15.  The  highest  order  mode 
is  determined  by  forcing  twice  the  higher  order  mode  spot  size  to  be  less  than  the 


30 


rod  radius,  2 WPti  <  Rr0d,  at  all  positions  inside  the  rod.  Once  the  mode  volume  is 
known,  then  the  number  of  ions  that  contribute  to  the  lasing  process  is  known,  and 
an  absolute  value  for  the  energy  and  power  can  be  determined. 

Finally,  to  find  the  spot  size  on  target  the  complex  beam  parameter  q  at  the 
outcoupler  needs  to  be  propagated  through  the  outcoupler,  and  then  propagated 
down  range  to  the  target.  Therefore,  the  transformation  matrix  through  the  outcou¬ 
pler  and  down  range  is  given  by 


A  B 
C  D 


'i 


0  1 


V  «M,  1 


(93) 


where  d  is  the  distance  to  the  target.  This  transformation  matrix  was  written  with 
the  assumption  that  the  index  of  refraction  of  the  outcoupler,  n=1.5,  the  outcoupler 
is  very  thin,  and  the  outside  surface  of  mirror  one  is  planar  (see  Figure  10).  Once 
the  spot  size  at  the  outcoupler  and  at  the  target  are  known,  then  the  divergence 
angle  can  be  found  by 


0Divcrg  =  arctan(W<orge^ — — )  (94) 
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2.5.3  Unstable  Cavity  For  the  stable  cavity,  it  was  possible  to  determine 
the  spot  size  everywhere  along  the  lasing  axis  and  thus  determine  mode  volume 
and  divergence  angle  of  the  beam.  For  the  unstable  cavity  it  is  very  difficult  to 
determine  the  spot  size  inside  the  cavity,  so  a  different  approach  will  be  used.  This 
approach  consists  of  two  steps.  The  first  step  will  be  to  find  the  cavity’s  round  trip 
magnification,  which  will  lead  to  finding  the  spot  size  at  the  outcoupler.  The  second 
step  will  be  to  equate  the  real  cavity  to  a  simplified  equivalent  cavity  to  determine 
the  divergence  angle  of  the  beam. 

The  round-trip  magnification  can  be  found  by  returning  to  equation  78.  This 
second-order  difference  equation  describes  the  rays  that  pass  through  any  cavity 
represented  by  an  ABCD  matrix.  Again,  the  solution  will  be  found  by  assuming  a 
solution  of  the  form  rj  =  r0FT.  Substituting  into  equation  78; 


raF 1 


F*-2^~^F+{AD-BC) 


=  0 


(95) 


In  order  for  there  to  be  a  non-trivial  solution,  the  expression  inside  of  the  brackets 
must  be  equal  to  zero,  and  therefore  F  can  be  found  by  using  the  quadratic  equation: 


F± 


A  +  D  , 
—± 


1 

2 


ABC- AD 


(96) 


Note  that  for  an  unstable  cavity  (^jp)2  >  AD  -  BC ,  and  so  the  solution  is  always 
real.  The  solution  F+  is  used  to  describe  a  positive- branch  unstable  cavity.  Positive 
branch  means  that  there  is  no  place  in  the  cavity  where  the  beam  comes  to  a  focus. 
The  solution  F. _  is  for  a  negative-branch  unstable  cavity.  Negative  branch  means 
that  there  is  one  or  more  internal  focal  point.  If  F+  and  F-  >  0,  then  the  cavity  is 
a  positive  branch,  and  if  F+  and  F-  <  0,  then  the  cavity  is  a  negative  branch.  The 
round-trip  magnification  for  an  unstable  cavity  is  given  by  M  =  and  therefore 
M  =  F.  The  fraction  of  photons  that  survive  after  one  round  trip  is  related  to  the 
magnification  by  T  =  and,  noting  that  F+  =  jr,  then  T+  =  Fl,  T_  =  and 
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the  fractional  unstable  out  coupling  S  =  1  -  T  (5). 

Now,  by  using  the  round  trip  magnification  and  the  diameter  of  the  outcoupler 
doc,  the  spot  size  at  the  outcoupler  can  be  found  by  woc  =  d0C\/~M. 

To  find  the  divergence  angle  of  the  beam  the  following  equations  are  used 
to  transform  the  real  cavity  into  a  simplified  equivalent  cavity  (10).  For  a  general 
single-ended  resonator  (full  round  trip): 

A  =  4fifi£f2  2  <7  2  1 

B  =  2g2L 

C  =  i{2gig2-9i  -gt) 

D  =  2g2  -  1  where 

91  =  1  “  ST 

92  =  1-% 

Knowing  ABCD  for  the  unstable  cavity  allows  the  determination  of  gi,  g2  and 
L  for  an  equivalent  unstable  cavity.  Then  from  (5) 
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(97) 


r  i 


[i  -  to)-1]*  -  i  +  gi1 

2  -  <7i_1  ~  92 1 


[1  ~  (gi^)"1]*  -  l  +  ft1 
2  -  gi1  ~  921 


(98) 


and  by  using  Figure  11, 


6  =  arctan 


doe(ri  + 1) 
2rir2£ 


(99) 
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III.  Validation  of  Model 


3.1  Analytic  Comparisons  to  Steady  State  Cases 

The  laser  model  can  be  broken  into  three  distinct  sections:  thermal,  rate  equa¬ 
tions,  and  cavity.  In  this  section  each  of  these  will  be  validated  by  comparing  them 
to  special  cases  which  have  analytic  solutions. 

The  first  section  to  be  validated  will  be  the  thermal  model.  By  taking  the  limit 
as  t  — >  oo  the  heat  model  will  go  to  a  steady  state  condition  ^  =  0,  and  equation  53 
can  be  written  as 

V2T  =  (100) 

fC 

d?T  ldT  _  q*(r) 

dr 2  r  dr  k 

allowing  q*(r)  =  Constant  —  q*  then 

1  (X)  =  -X 

dr  \  dr  )  k 

and  by  direct  integration 

T  =  -?£  +  Cln{r)  +  D  (103) 


Applying  the  boundary  conditions: 
T(0)  <  oo 

= few  - 

then  C  =  0  and 


^  =  i/t7mT-2TT-4fc- 

The  solution  can  be  written  as 


T(r)  =  T,am  +  £ 


R2-r2  + 


2  kR' 
h 


(104) 
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2  4  6  8: 

Time  (s) 

Figure  12.  Rod  Core  and  Wall  Temperature  Proceeding  to  the  Steady  State  Limit 

Figure  12  shows  the  temperature  distribution  of  the  rod  core  and  the  rod  wall 
as  a  function  of  time.  The  data  in  Figure  12  was  taken  by  using  a  flash  lamp  power 
of  10  kW,  a  power  fraction  for  heat  of  0.07,  a  pump  duration  of  10  s,  water  cooling 
with  a  velocity  of  10  m/s,  and  a  rod  length  of  0.1  m.  The  analytic  solution  for  the 
temperature  can  be  found  by  using  equations  61,  62,  and  104,  yielding  T(0)  =  436.90 
K  and  T(R)  =  394.05  K,  which  matches  the  laser  model  exactly. 

The  analytic  validation  for  the  rate  equations  will  proceed  in  a  similar  manner. 
As  t  — >  oo  then  =  0  and  equations  3  to  7  can  be  written  as 


$— (iV2  -  -NJj-  -$—  +  5  =  0 
n  Q\  lc  Trc 


(105) 


4-(iV2  -  ^JV,)/2  -  —  +  Hik  =  0 

n  g\  72i 


(106) 


crc.  ff2,ru  .  iV2/i  Ai 

$— {N2 - Ni)fi  +  — — - +  Rifi  =  0 

n  Q\  7*21  h  Tio 


(107) 
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Figure  13.  Continuous  Wave  Laser  Proceeding  to  Steady  State.  Shown  are  the 
many  relaxation  oscillations  before  steady  state  is  achieved. 
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and  substituting  equation  110  into  equation  109,  then  N?  can  be  written  as 


No  = 


hn _ 1_ 

lr  CTCTcni 


(112) 


Substituting  equation  112  into  equation  111  and  solving  for  $, 


AS  _  Df  _  lr  Tcav 
*  —  -^2^cav  » 


,  r  I - f- i?2rlo/lT' 

^2l/2  |_ 0"CTcav  *c 


(113) 


The  number  of  photons  outcoupled  into  the  beam  is  given  by 
and  has  units  of  photons  per  unit  volume  per  second,  and  the  outcoupled  power  is 
P  =  Using  equation  24  to  find  TcaV)  equation  49  to  find  r/av,  equation  15  to 

find  i?2>  equation  1  to  find  f\  and  /-j,  equation  41  to  find  <7,  and  also  using  Appendix 
C  for  T21,  no,  n,  v,  the  steady  state  power  per  unit  volume  can  be  found. 

As  an  example,  using  a  cavity  length  of  0.2  m,  a  mirror  1  reflectivity  of  0.9, 
a  mirror  2  reflectivity  of  0.98,  a  rod  length  of  0.1  m,  a  flash  lamp  power  of  10  kW 
and  a  power  fraction  for  lasing  of  0.03  (see  figure  13),  the  analytic  value  for  the 
power  is  83.4434  W/cc.  Comparing  that  to  the  numerical  value  = 

=  83.2649 W/cc  results  in  an  error  of  0.213  percent.  This  error  is  actually 
larger  than  the  real  error.  The  reason  is  that  the  final  output  of  the  peak  power  and 
mode  volume  from  the  model  are  truncated  which  produces  round-off  error.  If  this 
truncation  error  is  removed,  the  numerical  value  is  83.4473  W/cc,  and  the  error  is 
only  0.0047  percent. 

The  last  validation  will  be  for  the  cavity  by  examining  the  divergence  angle  as  it 
goes  to  the  far-field  limit.  Figure  14  is  a  plot  of  the  divergence  angle  for  a  stable  cavity 
versus  the  distance  from  the  outcoupler.  The  far-field  divergence  angle  is  given  by 
9  =  where  w0  is  the  radius  of  the  spot  at  the  beam  waist.  The  far-field  distance 
is  given  by  djj  >  y,  where- a  is  the  diameter  of  the  aperture,  and  A  is  the  wave 
length  of  the  laser  beam  (18).  Figure  14  was  generated  by  using  a  cavity  length  of 


38 


0.7 


0.6 
0.5 

Divergence 
0.4 
angle 

0.2 
0.1 
0 

0  5  10  15  20  25  30  35 

Distance  from  outcoupler  (m) 

Figure  14.  Divergence  Angle  as  it  Goes  to  the  Far  Field  Limit. 

0.2  m,  a  10  m  curvature  for  both  mirrors,  a  rod  length  of  0.1  m  and  by  varying  target 
distance  from  0  to  32  m.  Using  these  inputs,  the  laser  model  produced  a  spot  size  on 
the  outcoupler  of  0.5234  mm.  The  spot  size  at  the  beam  waist  w0  is  not  known,  but 
by  letting  wo  equal  the  spot  size  on  the  outcoupler,  a  lower  limit  on  the  divergence 
angle  can  be  found.  Calculating  the  divergence  angle  yields  9  =  0.6446  mrad.  The 
model  produces  a  divergence  angle  of  9  =  0.6519  mrad,  an  error  of  only  1  percent. 
To  calculate  the  far-field  distance,  the  diameter  of  the  outcoupling  aperture  must 
be  found.  Generally,  the  radius  of  the  outcoupling  aperture  must  be  twice  as  large 
as  the  spot  size  to  ensure  that  the  entire  transverse  mode  will  emerge.  Setting 
a  =  2.1296  mm,  djj  =  4.08  m,  which  appears  to  be  where  the  divergence  angle  in 
Figure  14  transitions  to  the  far-field  limit. 

3.2  Comparison  of  Literature  and  Modeled  Results 

With  a  realistic  model,  it  should  be  possible  to  match  experimental  data  fairly 
accurately.  The  problem  is  that  experimentalists  don’t  usually  give  enough  informa¬ 
tion  to  exactly  match  the  experiment.  This  requires  an  educated  ’’guess”  of  some  of 
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the  parameters. 

The  first  attempt  is  to  match  (11)  a  single  longitudinal  mode,  dye  Q-switched 
laser.  The  single  longitudinal  mode  was  selected  by  an  etalon  which  was  made  by 
placing  a  partially  transmissive  mirror  between  the  laser  rod  and  the  back  mirror. 
The  parameters  provided  by  the  article  are  as  follows: 

6.4mm  in  diameter  X  80mm  Nd:YAG  rod 

30J  flash  lamp  pulses 

100 fis  pump  duration 

Outcoupling  mirror  reflectivity  of  0.5 

Outcoupling  mirror  radius  of  curvature  6.7m 

Internal  mirror  reflectivity  0.85 

Internal  mirror  fiat 

Back  mirror  reflectivity  1 

Back  mirror  3  flat 

Cavity  length  of  77cm 

Etalon  length  of  25cm 

Total  length  (cavity  +  etalon)  of  102cm 

The  results  included  an  oscilloscope  trace  of  a  20ns  pulse  which  was  matched 
very  nicely  by  the  model  (see  figure  15).  Also  reported  was  a  10ns  pulse  with  an 
energy  of  7mJ.  By  changing  the  model’s  Q-switch  transmission  to  0.6  an  energy  of 
lOmJ  was  produced.  Therefore,  the  model  has  demonstrated  the  ability  to  predict 
single  mode,  Q-switched  pulse  shape  and  pulse  energy  fairly  accurately. 

The  next  attempt  will  be  to  match  the  results  of  a  repetitively  pulsed  laser. 
Experimental  data  from  (13)  was  reproduced  fairly  accurately  (See  figure  16). 

The  parameters  provided  by  the  article  are  as  follows: 

Water  cooling 

75mm  X  6mm  Nd:YAG  rod 

Cavity  length  of  600mm 
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Figure  15.  A  Single  Q-Switched  Pulse  that  Matches  an  Oscilloscope  Trace  from 
(11) 

Flat  mirrors 

2m  converging  lens  to  compensate  for  rod  lensing 

Outcoupling  mirror  reflectivity  of  0.60 

625Hz  repetition  rate 

160kW  input  power 

100^s  pump  duration 

The  reported  energy  in  the  first  pulse  was  35  mJ.  By  using  the  parameters 
provided,  and  by  varying  the  Q-switch  losses  and  the  laser  power  fraction  until  the 
first  pulse  was  approximately  35  mJ,  the  pulse  trend  V'as  matched  very  well.  The 
article  reported  a  12  percent  change  in  the  peak  power  between  the  1st  and  the  40th 
pulse,  and  the  model  produced  a  22  percent  change.  This  could  have  been  made  to 
match  the  experimental  results  more  closely  by  simply  decreasing  the  model’s  water 
flow  rate. 
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IV.  TVend  Analysis 


4-1  Energy  and  Power  vs.  Output  Coupling 

Figures  17  and  18  describe  the  peak  power  and  energy  of  a  Q-switched  laser 

versus  outcoupling,  for  various  input  power  settings.  Similar  plots  can  be  found  in 

(7)  and  they  help  to  demonstrate  that  the  model  is  working  properly.  Figures  17 

and  18  were  generated  by  varying  the  reflectivity  of  mirror  1  for  the  three  different 

* 

input  power  levels. 

4.2  High  Pulse  Repetition  Rate  Transient  Effects 

In  this  section  there  will  be  a  demonstration  of  the  model’s  ability  to  represent 
a  repetitively  pulsed  laser  by  exploring  trends  in  the  divergence  angle,  energy,  and 
peak  power. 

Figure  19  is  a  plot  of  the  divergence  angle  for  each  pulse  in  a  99  pulse  train, 
for  various  repetition  rates.  For  the  10  Hz  case,  there  is  a  large  increase  in  the 
divergence  angle,  and  then  it  levels  off.  The  reason  for  the  sudden  increase  is  that 
the  rod  is  heated  during  the  pump  cycle  and  with  the  extended  relaxation  period 
between  pulses  the  rod  has  a  chance  to  cool.  This  cooling  causes  a  radial  temperature 
gradient  in  the  laser  rod  which  is  responsible  for  producing  thermal  lensing.  The 
curve  flattens  out  because  the  rod  reaches  thermal  equilibrium,  i.e.,  the  temperature 
gradient  goes  to  a  constant,  and  therefore  the  lensing  of  the  rod  becomes  a  constant. 
The  100  Hz  case  proceeds  much  like  the  10  Hz  case  except  the  energy  is  deposited 
into  the  rod  at  a  higher  rate  which  causes  the  rod  to  heat  up  faster  than  the  10  Hz 
case.  This  increase  in  temperature  coupled  with  the  cooling  time  between  pulses 
allows  for  a  larger  thermal  gradient  to  develop.  Finally,  the  1000  Hz  case  initially 
rises  much  slower  than  the  two  previous  cases.  This  is  because  a  large  number  of 
pulses  occur  before  a  significant  amount  of  heat  can  be  removed  from  the  rod  wall 
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Figure  17.  Energy  Out  vs.  Outcoupling  for  a  Single  Pulsed  Q-Switched  Laser  at 
Various  Input  Power  Settings 
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Figure  18.  Peak  Power  vs.  Outcoupling  for  a  Single  Pulsed  Q-Switched  Laser  at 
Various  Input  Power  Settings 


44 


0  10  20  30  40  50  60  70  80  90  100 

Laser  Pulse 

Figure  19.  Divergence  Angle  vs.  Laser  Pulse  for  10  J  Input  Energy  and  Various 
Repetition  Rates. 

and  thus  initially  a  significant  thermal  gradient  is  not  produced.  Note,  the  average 
temperature  of  the  laser  rod  is  higher  for  the  1000  Hz  case  than  it  is  for  the  other 
two  cases.  This  is  because  the  total  amount  of  energy  deposited  into  the  rod  is  the 
same  for  all  cases,  but  in  the  10  Hz  and  the  100  Hz  cases,  there  is  more  time  to  cool 
between  pulses  than  for  the  1000  Hz  case. 

Figure  19  was  generated  by  using  a  pump  duration  of  100  fj,s  and  a  flash  lamp 
power  of  100  kW.  The  relaxation  time  was  set  to  0  [is  for  the  100  Hz  case,  to  9900 
fis  for  the  100  Hz  case,  and  to  99900  fis  for  the  10  Hz  case. 

Figure  20  very  graphically  shows  the  effect  of  cooling  the  laser  rod.  The  top 
curve  is  the  first  pulse  in  a  40  pulse  train;  it  is  the  same  for  all  three  cooling  rates. 
Noting  the  trend  in  Figure  20,  it  can  be  seen  that  as  the  coolant  flow  rate  is  increased, 
the  power  in  the  40th  pulse  decreases.  This  means  that  the  higher  the  flow  rate,  the 
faster  heat  is  removed  from  the  surface  of  the  laser  rod,  thus  causing  a  larger  radial 
thermal  gradient.  This  causes  lensing  in  the  rod  which  serves  to  decrease  the  mode 
volume  and  therefore  reduces  the  number  of  ions  that  contribute  to  the  laser  beam. 
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Figure  20.  Trend  of  i'.iQ  1st  Pulse  and  the  40th  Pulse  for  Three  Different  Coolant 
Velocity  Settings.  The  first  pulse  is  the  same  for  all  fluid  velocities. 

The  total  energy  into  the  laser  rod  is  a  constant,  and  for  the  longer  pulses,  the  flash 
lamp  power  decreases  and  therefore  the  pump  rate  decreases.  This  lower  pump  rate 
results  in  a  smaller  inversion  even  though  the  flash  lamp  is  on  for  a  longer  time. 

Figure  20  was  made  by  varying  pump  duration,  Q-switch  delay,  relaxation  time 
and  flash  lamp  power  to  keep  the  total  input  energy  constant  at  10  J  and  to  hold 
the  repetition  rate  constant  at  1000  Hz. 

Figure  21  is  a  plot  of  the  peak  power  of  the  40th  pulse  versus  repetition  rate. 
It  was  generated  by  using  a  pump  duration  of  1000  /is,  a  Q-switch  delay  of  999.5  /is, 
and  by  varying  the  relaxation  time  from  0  to  9000  /is  for  15,  25,  and,  35  kW  input 
powers. 

What  can  be  seen  from  Figure  21  is  that  for  very  low  repetition  rates  the 
laser  rod  has  a  chance  to  fully  cool  between  laser  pulses,  but  as  the  repetition  rate 
increases  there  is  less  time  for  the  laser  rod  to  completely  cool,  allowing  a  thermal 
gradient  to  develop  inside  the  rod.  This  thermal  gradient  is  responsible  for  producing 
lensing  which  causes  a  reduction  in  the  mode  volume  and  thus  a  reduction  in  the 
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Figure  21.  Peak  Power  vs.  Repetition  Rate  for  Various  Input  Powers  and  Constant 
Pump  Duration  of  1  ms 

peak  power.  For  higher  repetition  rates,  the  40th  pulse  is  produced  before  the  laser 
rod  develops  a  large  thermal  gradient,  which  translates  into  only  a  small  decrease 
in  mode  volume.  For  the  lower  power  settings  the  peak  power  is  less  sensitive  to 
the  repetition  rate  because  there  is  less  heating  of  the  rod  and  only  a  small  thermal 
gradient  can  develop. 

The  next  case  will  examine  the  energy  in  the  40th  pulse  of  a  continuously 
pumped,  Q-switched  laser  (see  Figure  22).  At  the  low  repetition  rates,  there  is  a 
significant  amount  of  thermal  lensing,  because  the  lamp  has  been  on  for  a  long  period 
of  time  allowing  a  large  thermal  gradient  to  develop  which  causes  a  reduction  in  mode 
volume.  As  the  pulse  rate  increases,  the  40tli  pulse  can  emerge  from  the  cavity  before 
the  rod  can  develop  a  significant  thermal  gradient.  Once  again,  the  curve  begins  to 
fall  off  because  there  isn’t  enough  time  between  pulses  to  fully  populate  the  upper 
laser  level,  causing  a  decrease  in  the  output  energy.  The  sudden  rise  is  from  inversion 
left  over  from  the  39th  pulse.  If  the  inversion  isn’t  high  enough  for  the  39th  pulse  to 
lase  then  the  inversion  will  contribute  to  the  40th  pulse.  The  sudden  drop  is  because 
the  40th  pulse  has  become  the  pulse  without  enough  inversion  to  lase. 
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Figure  22.  Energy  in  the  40th  Pulse  vs.  Repetition  Rate. 

Figure  22  was  generated  by  using  a  flash  lamp  power  of  10  kW,  a  relaxation 
time  of  0,  and  by  varying  the  pump  duration  from  100  to  10000  ns.  A  relaxation 
time  of  0  means  that  the  flash  lamp  will  not  be  turned  off  between  pulses. 


V.  Conclusion 


5.1  Conclusion 

By  modeling  the  laser  rate  equations,  the  cavity  design,  and  the  temperature 
distribution  in  the  rod,  it  is  possible  to  do  an  accurate  job  of  describing  an  Nd:YAG 
laser  for  both  continuous  wave  and  pulsed  operation.  Each  piece  is  necessary  to 
accurately  represent  the  laser.  The  rate  equations  are  responsible  for  producing  the 
detailed  laser  pulse  shape  and  for  producing  the  outputs  of  power  and  energy.  The 
rate  equations  are  limited  in  that  they  can  only  produce  energy  and  power  on  a 
per  unit  volume  basis.  The  cavity  model  is  included  to  calculate  the  fraction  of  rod 
volume  that  contributes  to  the  beam,  making  absolute  energy  and  power  calculations 
possible.  To  model  the  cavity  properly,  the  thermal  lensing  of  the  laser  rod  needed  to 
be  included.  This  lensing  has  a  significant  effect  on  the  mode  volume  and  therefore 
has  a  significant  effect  on  the  energy  and  power  predicted  by  the  model. 

5.2  Suggested  Improvements 

While  coding  the  laser  model,  a  QuickBASIC  limit  was  discovered.  QuickBA¬ 
SIC  is  limited  to  a  maximum  module  size  of  64K,  which  ultimately  prevented  the 
addition  of  new  material  to  the  laser  model.  In  order  to  allow  future  expansion  of 
the  model,  the  code  would  need  to  be  broken  into  smaller  modules. 

After  the  code  is  broken  into  smaller  modules,  the  first  addition  should  be  a 
more  realistic  Q-switch.  Presently  the  Q-switch  goes  from  off  to  on  instantaneously 
and  stays  open  until  the  end  of  the  integration  of  the  pulse.  A  more  realistic  Q-switch 
should  include  a  finite  opening  and  closing  time. 

The  present  unstable  cavity  code  does  not  give  any  information  about  the  beam 
inside  the  cavity.  This  lack  of  information  means  that  an  accurate  evaluation  of  the 
mode  volume  can  not  be  determined.  Presently  a  mode  volume  equal  to  the  size  of 
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the  laser  rod  is  assumed.  This  is  not  a  good  assumption  and  therefore  the  energy 
and  power  for  the  unstable  cavity  design  is  somewhat  questionable. 

A  very  difficult  addition  would  be  a  realistic  flash  lamp  model.  This  addition 
is  not  needed  for  the  Nd  ion  laser  but,  if  there  is  a  future  desire  to  extend  the  code 
to  other  laser  materials  then  this  would  be  an  important  modification. 

The  last  set  of  improvements  are  all  related  to  temperature  dependence.  The 
first  is  thermally  induced  birefringence.  As  the  rod  develops  a  radial  temperature 
gradient,  stress  is  produced  inside  the  laser  rod.  This  stress  produces  birefringence 
effects.  Thermally  induced  birefringence  is  as  much  a  design  consideration  as  the 
thermal  lensing  of  the  laser  rod  and  should  be  added  to  more  completely  model  the 
Nd  laser. 

,  The  stress  that  produces  the  birefringence,  in  extreme  cases,  can  also  cause 
the  laser  rod  to  fracture.  These  laser  rod  stress  limits  should  be  added  to  make 
very  high  energy  laser  pulse  evaluations  more  realistic.  Another  thermal  effect  is  the 
laser  rod  expansion  along  the  lasing  axis.  This  lengthening  of  the  rod  affects  the  laser 
cavity’s  design  and  therefore  the  cavity’s  stability.  And  finally,  temperature  changes 
in  the  laser  rod  produce  shifts  in  the  laser  frequency.  Tracking  the  changes  in  the 
frequency  would  be  useful  in  matching  the  laser  output  to  detectors  in  applications 
such  as  laser  radar. 
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Appendix  A.  Numerical  Techniques 


The  model  uses  an  Adams  predictor-corrector  method  to  solve  the  differential 
equations  in  chapter  2.  Because  this  method  requires  a  set  of  initial  data  points,  a 
fourth  order  Runge-Kutta  integration  routine  is  used  to  start  the  Adams  predictor 
corrector  routine. 

The  fourth  order  Runge-Kutta  technique  uses  the  following  algorithm  (3): 

Define  the  rate  equations  as  functions  to  be  called  by  the  integration  routines. 

Fi{<f> ,  m,  n2,  n3)  =  ^(n2  -  — «i)  -  (f>  +  s  (114) 

9i 


F2(^,rii,n2,n3)  =  <f>(n2  -  —nx)fXy  +  y-j-  +  —  fx  -  y-  +  rxfx 

9 1  ‘r  ‘21 J  2  ‘31  ‘10 


(115) 


jRiW,ni,n2)n3)  =  -0(n2  -  — n,)/2|l  +  ™/2 

9 1  *r  ‘32 


n2 

<2 


+  r2/2 


(116) 


F4(<£,ni,n2,n3) 


n3  . 

— +  ^3 
‘3 


(117) 


Starting  with  the  following  initial  conditions 


0  <t<b  m(0)  =  0  n2(0)  =  0  rc3(0)  =  0  ^(0)  =  0  (118) 

where  b  =  the  time  to  integrate  to  and  St  =  the  incremental  time  step. 

Set  N  =  i 
Set  t  =  0 
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Loop  N  times  to  integrate  to  time  b 

Loop  i  from  1  to  4 
Set  K\{  -  ^jPi(^,ni,n2,n3) 

Loop  i  from  1  to  4 

Set  It2i  —  6tFi(<f>+  ^-,711  +  ^,n2  +  ^-,713  +  ^-) 

Loop  i  from  1  to  4 

Set  Kzi  —  StFi((f>  +  ^-,711  +  ^-^2  +  ^ip,7i3  +  ^-) 

Loop  i  from  1  to  4 

Set  I<4{  =  6tFi((j)  +  Kzi )  77 1  +  Kz2i  n2  +  It 33>  72.3  -f  K34) 

Set  </)'-({>+  Zi+M&Mtlk 

Set  7i j  =  7ix  +  {h±M2±MtK 1 
Set  n'2  =  n2  -f 

Set  n'3  =  m  +  &±2-&±!ia)+K* 
set  t’  =  t  +  St 

The  primed  quantities  are  the  values  at  the  time  t\ 

The  Runge-Kutta  routine  is  very  capable  and  could  have  been  used  exclusively 
to  integrate  the  set  of  differential  equations.  The  problem  with  the  routine  is  that  it 
is  very  slow.  Also,  it  is  only  a  fourth-order  routine,  which  means  that  a  very  small 
time  step  is  needed  to  accurately  integrate  the  rate  equations.  The  small  time  step 
translates  into  a  very  large  number  of  points,  a  lot  of  computation  time  to  calculate 
the  points,  and  a  lot  of  storage  space  to  store  the  points.  To  get  around  these 
limitations,  the  Runge-Kutta  routine  is  only  used  to  initiate  a  sixth-order  Adams 
predictor-corrector  routine.  The  sixth-order  routine  allows  a  larger  time  step  to  be 
used  and  still  achieve  the  same  accuracy  as  the  Runge-Kutta  routine.  Also,  the 
Adams  method  is,  by  design,  a  faster  technique.  Therefore,  there  is  a  two-fold  speed 
increase  by  employing  this  routine. 
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There  is  a  problem  with  initiating  a  sixth-order  routine  with  a  fourth-order 
routine.  The  lower  accuracy  of  the  starting  routine  could  introduce  error  into  the 
integration.  This  problem  is  overcome  by  running  the  Runge-Kutta  routine  with  a 
time  step  of  j  and  using  every  other  point  generated  to  start  the  Adams  routine. 

The  following  is  the  Adams  predictor-corrector  algorithm  (4): 

Call  the  Runge-Kutta  routine  10  times  with  a  time  step  of  j  and  use  the  0th,  2nd, 
4th,  6th,  8th,  and  10th  points  to  start  the  Adams  routine. 

Set  N  =  i 

Set  t  =  5 St 

Loop  i  from  6  to  N  to  integrate  to  time  b 
Predict 

Set  4/  =  4m  +  -  gi :F,(i  ■ -  St)  +  H  Ji(‘  -2 st)~  HF,(t  -  3«)  + 

Set  »;  =  n,m  +  -  g|F2((  -  St)  +  ?§!§F2(i  -  2 St)  -  »F2(i  - 

3«)  +  -  4St)  -  ISsW*  -  5ft)] 

Set  n'2  =  n2(i)  +  ft|^F»(t)  -  Zf§F3(t  -  St)  +  «§F,(i  -  2 St)  -  gf  F3(t  - 
3«)  +  ggFs(<  ~  ~  flkW  ~  5«)J 

Set  n't  =  «,(*)  +  «[®R,(t)  -  -  St)  +  ffgF4(«  -  2  St)  -  f «F((  - 

3«)  +  g gF.(<  "  4 St)  ~  F„(«  -  5ft)] 

Correct 

Set  4"  =  m  +  ft[£|fi(t  +  St)  +  USf.fi)  -  &F,(«  -  «)  +  fskW  - 2{<)  - 
SJ?.(‘-3«)  +  I5of.(*-4«)] 

Set  <  =  n,(t)  +  Stl£LF,(i +«)  +  ^F,(i)  -  g|F2(t  -  ft)  + - 2 St) - 

S«-»)  +  S  F»(«-4«)] 

Set  n'i  =  n2(i)  +  «[£&F,(<  +  St)  +  Jg=F,(()  -  ^ F,(t  -  St)  +  F3(t  -  2 St)  - 

Jg^,(i-3«)+TgiF,(i-4«)] 
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SetnJ  =  n3(t)+ft[^(‘  +  «)  +  ^i^(i)-S^(<-«)  +  S^(i-2«)- 
l£0F4(t  -  3 St)  +  ^F,(t  -  4«)1 

set  t’  =  t  +  6t 

The  primed  quantities  are  the  predicted  values  and  the  double  primed  quanti¬ 
ties  are  the  corrected  values  at  time  t\  A  useful  benefit  from  the  routine  is  that  the 
corrector  portion  is  closed-form,  making  it  possible  to  iterate  the  corrector  section  to 
increase  the  accuracy.  The  laser  model  utilizes  this  feature  to  produce  better  results. 
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Appendix  B.  Code  Instructions 


B.l  User  Input 

Insert  floppy  into  drive  A: 
type  A:  return 
type  Nd  return 

Select  the  number  pad  operational 

Push  2  to  move  the  highlight  bar  up  and  push  8  to  move  the  highlight  bar 

down. 

Press  return  to  select  the  highlighted  option. 

Again  press  2  and  8  to  move  up  and  down  until  the  proper  choice  is  highlighted 
and  press  return.  Input  the  new  value  at  the  prompt  and  proceed  as  above  to  the 
next  value  to  change. 

When  all  parameters  are  set  then  highlight  RUN  MODEL  and  press  return  to 
run  the  code. 

CAVITY 

Cavity  Length 

Must  be  longer  than  the  laser  rod 
Mirror  1 

This  is  the  reflectivity  of  the  outcoupling  mirror.  The  reflectivity  is  restricted  to  the 
range  of  0.1  to  0.99999. 

Mirror  2 

This  is  the  reflectivity  of  the  back  mirror  and  is  restricted  to  the  range  of  0.1  to 
0.99999. 

Q-Switch  Trans 

This  is  the  single-pass  transmission  of  the  Q-Switch  and  ranges  from  0.1  to  0.99999. 
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Curvature  Mirror  1 

The  radius  of  curvature  of  the  outcoupling  mirror  and  its  value  can  range  from  -1E10 
to  +1E10. 

Curvature  Mirror  2 

The  radius  of  curvature  of  the  back  mirror  and  its  value  can  range  from  -1E10  to 
+1E10. 

Diameter  Mirror  1 

Diameter  of  the  outcoupling  mirror.  Only  used  during  unstable  cavity  calculations. 
Its  value  must  be  greater  than  zero. 

Target  Distance 

This  is  the  distance  to  some  arbitrary  target  and  is  used  to  calculate  the  laser  spot 
size  on  target. 

ROD-COOLANT 

Rod  Length 

Must  be  less  than  the  length  of  the  cavity  and  greater  than  zero. 

Rod  Radius 

The  radius  of  the  laser  rod,  and  must  be  greater  than  zero. 

Curvature  Rod  Ends 

The  radius  of  curvature  at  the  ends  of  the  laser  rod.  Can  range  from  -1E10  to  +  1E10. 
Laser  Material 

This  is  the  material  the  laser  rod  is  made  of  where  1  =  Nd:YAG,  2  =  Nd:Silicate 
”ED-2”,  and  3  =  Nd:Phosphate  ”Q-88”. 

Percent  Nd:  Ions 

Percentage  of  Nd:  Ions  in  the  laser  rod  and  must  be  greater  than  zero. 

Coolant  Velocity 

The  velocity  that  the  cooling  fluid  flows  over  the  laser  rod.  Its  value  must  be  greater 
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Figure  23.  Flash  Lamp  Duty  Cycle 

than  zero. 

Coolant 

This  selects  the  cooling  fluid  where  1  =  Water  and  2  =  Dry  Air. 

Coolant  Temperature 

This  is  the  bulk  temperature  of  the  cooling  fluid.  It  is  also  the  temperature  of  the 
laser  rod.  The  coolant  temperature  range  is  different  for  the  two  different  cooling 
fluids.  For  water  the  temperature  range  is  from  273  K  to  600  K,  and  for  air  the 
temperature  range  is  from  223  K  to  2500  K. 

FLASH  LAMP 

Pump  Duration 

The  length  of  time  that  the  flash  lamp  is  turned  on  (see  Figure  23).  NOTE:  (Pump 
Duration  +  Relaxation  Time)  must  be  greater  than  or  equal  to  the  (Q-Switch  Delay 
+  Integrate  to). 

Q-Switch  Delay 

The  delay  from  the  time  the  flash  lamp  is  turned  on  to  the  time  the  Q-switch  is 
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opened.  By  setting  this  value  to  zero  the  losses  from  the  Q-switch  are  neglected. 
This  has  the  same  effect  as  removing  the  Q-switch  from  the  cavity.  NOTE:  (Q-Switch 
Delay  +  Integrate  to)  must  be  less  than  or  equal  to  (Pump  Duration  +  Relaxation 
Time). 

Relaxation  Time 

The  length  of  time  from  flash  lamp  turn-off  to  flash  lamp  turn-on  (see  Figure  23). 
NOTE:  (Pump  Duration  +  Relaxation  Time)  must  be  greater  than  or  equal  to  the 
(Q-Switch  Delay  +  Integrate  to). 

Flash  Lamp  Power 

The  electrical  power  into  the  flash  lamp  or  diode  pump  array.  This  value  must  be 
greater  than  zero. 

Power  Fraction  Laser 

The  fraction  of  flash  lamp  electrical  power  that  goes  into  the  laser.  A  typical  value 
for  a  flash  lamp  is  about  0.03,  and  a  typical  value  for  a  diode  pump  source  is  about 
0.2. 

Power  Fraction  Heat 

The  fraction  of  flash  lamp  electrical  power  that  goes  into  heating  the  laser  rod.  A 
typical  value  for  a  flash  lamp  is  about  0.07,  and  a  typical  value  for  a  diode  pump 
source  is  about  0.2. 

INTEGRATION 

Number  of  pulses 

The  number  of  laser  pulses  to  be  integrated.  The  value  must  be  an  integer  ranging 
from  1  to  100. 

Integrate  to 

” Integrate  to”  always  starts  from  the  time  the  Q-Switch  is  opened  and  is  the  amount 
of  time  the  integrator  will  integrate.  This  value  must  be  greater  than  zero  and  is 
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typically  set  between  .5  and  1  for  most  Q-Switched  pulses,  but  can  range  as  high  as 
100  to  150  for  CW  operation. 

Step  Size 

This  is  the  incremental  step  s'ze  that  the  integrator  will  use  to  integrate  the  laser 
pulse.  The  smaller  the  step  size,  the  more  accurately  the  laser  pulse  is  integrated, 
but  the  cost  is  a  longer  wait  for  the  output.  Also,  the  smaller  the  step  size,  the  larger 
the  number  of  points  needed  to  integrate  the  pulse.  This  value  should  not  be  greater 
than  0.7  for  CW  operation  and  not  greater  than  0.2  for  Q-Switch  operation.  These 
limits  will  help  eliminate  the  OVERFLOW  error  by  the  computer. 

File  Storage  Path 

This  is  the  path  name  the  model  will  use  to  locate/save  the  four  data  files  created, 
i.e.  ”C :/”,  ”D:/”,”C:/SUBDIRECTORY/”  etc.  Disregard  the  quotation  marks  in 
these  examples. 

Mode  of  Operation 

This  is  a  simple  switch  that  determines  Q-Switch  operation  versus  Continuous  Wave 
operation.  1  =  Q-Switch  and  2  =  CW.  NOTE:  When  CW  is  selected,  the  heat 
model  just  produces  the  steady  state  temperature  distribution  inside  the  laser  rod. 
Also,  when  CW  is  selected  (Pump  Pulse  =  Integrate  to),  the  Q-switch  delay  and  the 
relaxation  time  are  set  to  0. 

Store  Power  File 

This  is  a  simple  switch  that  determines  if  the  power  file  will  be  saved  to  disk;  1  = 
Yes  and  2  =  No.  When  a  large  number  of  pulses  are  integrated,  the  power  file  can 
become  quite  large,  taking  up  a  lot  of  storage  space,  and  storing  the  power  file  is 
not  advisable.  If  the  power  file  is  not  saved  then  the  two  power  plots  can  not  be 
displayed. 

Number  Of  Points 

This  is  not  an  input  but  by  setting  the  (Step  Size)  and  (Integrate  to)  the  value 
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changes.  This  value  should  not  exceed  10000  for  a  single  pulse.  Exceeding  10000 
points  in  a  single  pulse  will  produce  a  SUBSCRIPT  OUT  OF  RANGE  error. 

B.2  Model  Output 

After  running  the  model  the  default  output  is  information  about  each  laser 
pulse.  To  continue  on  past  this  point  press  Return  until  the  main  menu  returns. 
Once  the  model  has  been  run,  a  new  menu  item  appears,  PLOT  DATA.  By  selecting 
PLOT  DATA  another  menu  appears.  The  first  entry  is 

Temp  vs  Time 

This  will  produce  a  plot  of  the  rod  core  and  the  rod  wall  temperature  with  respect 
to  time.  To  exit  back  to  the  menu,  press  Return  twice. 

Power  vs  Time 

This  is  a  plot  of  laser  power  with  respect  to  time  for  all  laser  pulses.  To  exit  back 
to  the  menu  pre^s  Return  twice.  This  will  not  be  displayed  if  the  store  power  file 
option  was  set  to  NO. 

Pulse  P  vs  T 

This  is  a  plot  of  the  laser  power  for  each  individual  pulse.  To  exit  back  to  the  menu 
press  Return  once  and  answer  no.  This  will  not  be  displayed  if  the  store  power  file 
option  was  set  to  NO. 

Pulse  Data 

This  is  the  same  information  that  is  displayed  when  the  code  finishes  running  a  case. 
Below  is  a  description  of  what  these  outputs  mean. 

Shell  to  Dos 

This  option  allows  the  user  to  temporarily  exit  the  code  and  use  DOS  commands. 
When  the  user  is  finished  the  model  can  be  reentered  by  typing  exit  followed  by 
return. 
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The  following  is  a  description  of  the  information  displayed  by  the  Pulse  Data 
option. 

PULSE 

Identifies  the  pulse  number. 

STABILITY 

This  has  three  possible  outputs:  Stable,  Positive,  and  Negative.  Stable  means  that 
the  laser  cavity  was  stable  for  the  laser  pulse,  Positive  means  positive  branch  unstable 
cavity,  and  Negative  means  negative  branch  unstable  cavity. 

STABILITY 

This  is  the  numerical  value  for  stability  and  generally  must  be  less  than  two  for  the 
cavity  to  be  stable. 

ENERGY  0,0 

The  energy  in  the  laser  pulse  for  the  0,0  transverse  mode.  The  units  are  Joules. 
PEAK  PWR  0,0 

The  peak  power  of  the  laser  pulse  for  the  0,0  transverse  mode.  The  units  are  Watts. 
ENERGY  p,0 

The  energy  in  the  laser  pulse  for  the  p,0  transverse  mode.  The  units  are  Joules. 
PEAK  PWR  p,0 

The  peak  power  of  the  laser  pulse  for  the  p,0  transverse  mode.  The  units  are  Watts. 
FWHM 

Full-width  at  half-maximum  of  the  laser  pulse  power.  The  units  are  seconds. 
PULSE  ST/  RT 

This  is  the  real-ti;  'tart  of  the  laser  pulse,  and  is  used  to  help  isolate  the  pulse 
when  plotting.  This  value  is  meaningless  for  CW  operation. 

PULSE  END 

This  is  the  real-time  end  of  the  laser  pulse,  and  is  used  to  help  isolate  the  pulse  when 
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plotting.  This  value  is  meaningless  for  CW  operation. 

TEMP  r=0 

The  temperature  of  the  laser  rod  core  during  the  laser  pulse  in  Kelvins. 

TEMP  r=0 

The  temperature  of  the  laser  rod  wall  during  the  laser  pulse  in  Kelvins. 

SPOT  Ml 

Spot  size  radius  on  the  outcoupling  mirror  in  millimeters. 

SPOT  M2 

Spot  size  radius  on  the  back  mirror  in  millimeters. 

SPOT  min 

Minimum  spot  size  radius  inside  the  laser  rod  in  millimeters. 

SPOT  max 

Maximum  spot  size  radius  inside  the  laser  rod  in  millimeters. 

MODEV  0,0 

Volume  of  the  laser  rod  used  by  the  0,0  laser  mode  in  cubic  centimeters. 

MODE  p,0 

Indicates  the  highest  order  transverse  mode  that  the  laser  can  operate  in  where 
p=mode. 

MODEV  P,0 

Volume  of  the  laser  rod  used  by  the  P,0  laser  mode  in  cubic  centimeters. 

SPOT  oc 

Spot  size  radius  at  the  outcoupling  mirror  outside  the  cavity  in  milbrneters. 

SPOT  TARGET 

Spot  size  radius  on  target  in  meters. 

SSGain 

Small-signal  gain  in  the  laser  rod  at  the  start  of  the  Q-switched  pulse  in  cm'1. 
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UNSTAB  OUT 

Fractional  outcoupling  for  an  unstable  cavity. 

Dvrg  anle 

Beam  divergence  angle  in  milliradians. 

Rod  Condition 

There  two  possible  outputs  including  GOOD  and  MELTED. 

There  are  also  four  output  files  stored  on  disk.  These  files  are  stored  according 
to  the  path  specified  on  the  input  screen.  The  first  file  is  titled  ” Input.dat”,  this 
file  contains  a  snapshot  of  the  input  screen.  The  second  file  is  titled  ” Pulse.dat” 
and  contains  the  same  information  displayed  by  the  pulse  data  option.  The  third 
file  is  named  ’’TmpTme.dat”  and  contains  the  laser  rod  temperature  information. 
The  information  is  stored  in  the  format  Time  (ms),  Rod  Core  Temp  (K),  Rod  Wall 
Temp  (K).  The  last  file  is  named  ’’Pow.dat”  and  holds  the  power  information  in  the 
format  Time  (ms),  Power  (kW). 


63 


Appendix  C.  Table  of  Constants 


Table  6.  Laser  Material  Constants  along  with  the  Source  Reference. 


Materieal 

Nd:YAG 

NdrSil  ED-2 

NdtPhos  Q-88 

Scatter  losses  by  rod  (/m) 

.2 

.3 

.3 

Einstein  A21  coefficient  (/s) 

820 

1380 

1941  1 

T10  (s) 

HKmfiflHi 

MfrlfH'Blfrl 

UBBII 

T20  (s) 

00 1 

00 1 

00 1 

T21  (s) 

mttmiam 

kszssih 

T30  (s) 

00  1 

00 1 

00  1 

T31  (s) 

00  1 

00 1 

00  1 

T32  (s) 

■EHSilHI 

NHHSHiJ 

Index  of  Refraction 

1.82 

1.5554 

1.5310 

Laser  Frequency  Hz 

^mmsusM 

3 

1 

s 

Line  Width  Hz 

1.4985Z1011 

Energy  Levis  in  Manafold  1 

6 

2 

2 

Energy  of  level  1  /cm 

0 

0 

0 

Degeneracy  of  level  1 

2 

4 1 

4 1 

Energy  of  level  2  /cm 

26 

310  1 

Degeneracy  of  level  2 

2 

8  1 

8  1 

Energy  of  level  3  /cm 

108 

_ 

- 

Degeneracy  of  level  3 

2 

- 

- 

Energy  of  level  4  /cm 

146 

- 

- 

Degeneracy  of  level  4 

2 

- 

- 

Energy  of  level  5  /cm 

459 

- 

- 

Degeneracy  of  level  5 

2 

- 

Energy  of  level  6  /cm 

512 

- 

- 

Degeneracy  of  level  6 

2 

- 

- 

Energy  level  transitioned  to 

3 

1 

1 

1Warning:  An  experimental  value  could  not  be  found.  The  value  given  is  only  a  ’’reasonable” 
approximation  of  the  real  value. 
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Table  7.  Laser  Material  Constants  alone  with  the  Source  Reference  Continued. 


Materieal 

Nd:YAG 

Nd:Sil  ED-2 

Energy  Levis  in  Manafold  2 

2 

2 

Energy  of  level  1  /cm 

0 

0 

Degeneracy  of  level  1 

2 

2 

Energy  of  level  2  /cm 

84 

110  1 

Degeneracy  of  level  2 

2 

2 

Energy  level  transitioned  from 

2 

1 

Conductivity 

13 

1.255 

590 

920.48 

Density 

4560 

2510 

Melting  Point  K 


Source 


2243 


(5)  (7)  (8) 
(14)  (15) 


1000  1 


(7)  (16) 
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